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Models consisting of chains of particles that are 
coupled to their neighbours appear in many applica- 
tions in physics or engineering, such as in the 
study of dynamics of mono-atomic and multi-atomic 
lattices, the resonances of crystals with impurities 
and the response of damaged bladed discs. Analytical 
properties of the dynamic responses of such disturbed 
chains of identical springs and masses are presented, 
including when damping is present. Several remark- 
able properties in the location of the resonances 
(poles) and anti-resonances (zeros) of the dis- 
placements in the frequency domain are presented 
and proved. In particular, it is shown that there 
exists an elliptical region in the frequency-disturbance 
magnitude plane from which zeros are excluded and 
the discrete values of the frequency and disturbance at 
which double poles occur are identified. A particular 
focus is on a local disturbance, such as when a spring 
or damper is modified at or between the first and 
last masses. It is demonstrated how, notably through 
normalization, the techniques and results of the paper 
apply to a broad category of more complex systems in 
physics, chemistry and engineering. 



1. Introduction 

The dynamics of collinear chains of particles has 
been studied for a long time. Although apparently 
simple, these systems exhibit remarkable properties. 
Their interest is not limited to the academic realm as they 
model the fundamental vibrations of crystals in solid- 
state physics [1,2], the atomic and molecular dynamics 
of chains of molecules in physics, chemistry and biology 
[3-7], as well as the behaviour of real-life objects such 
as structures with repetitive components or rods and 
beams that are widely used in engineering. They are 
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directly related to linear algebra in which they are represented by particular matrices, such as 
tri-diagonal, multi-diagonal, circulant or Toeplitz matrices, as well as to the theory of orthogonal 
polynomials. They are also often used to study or illustrate theoretical aspects like the effect of 
disorder on systems [8,9]. Their properties indicate how the systems can behave when used within 
a controlled system or when some of their components are perturbed or damaged, for example in 
the case of damaged bladed discs or impurities in crystals [10-15]. 

However simple these systems appear, however important they are and for however long they 
have been studied, they are still under investigation. The present contribution is notably closely 
related to recent work on analytical expressions of the eigenvalues, eigenvectors and inverse of tri- 
diagonal matrices that have two or four of their corner coefficients disturbed [16-21]. The focus 
here is on the location of not only poles but also zeros of the transfer functions of a disturbed 
system in the frequency-disturbance magnitude plane. These loci are important system properties 
as they notably indicate if a system is stable, observable or controllable [22-24]. The loci of zeros 
also allow the localization of stable and unbiased frequency points in the theory of uncertain 
dynamic systems [13,25]. Analytical explicit expressions of the transfer functions of the system 
are derived by making use of Chebyshev polynomials and low-rank updates. 

After presentation of the system and resulting equations of motion in §2, the document is 
mainly split into the study of the nominal or undisturbed system, and that of the system with 
rank-one disturbance. For both, the exact expressions of all transfer functions are first presented 
in §3. Properties of the nominal system are then presented in §4. This includes the explicit 
expressions for the locations of poles and zeros of the transfer functions. Similar results and 
properties are presented in §5 for the disturbed system. Along the way, remarkable properties 
of the transfer function from the first to last mass are presented. In the case of the nominal 
system, it is shown in §4c that, besides the known property that this transfer function has no 
zero, its magnitude is also excluded from a circular region. Similarly, for the disturbed system, 
it is shown in §5b(i) that no zeros of this transfer function exist in an elliptical region of the real 
plane defined by the frequency parameter and the disturbance magnitude. Other exclusion and 
inclusion properties for the loci of poles and zeros are presented, such as the location of multiple 
poles of the disturbed system in §5#(ii). 

2. Nominal and disturbed spring-mass systems 

Considered is a nominal collinear system of alternating springs and masses, as illustrated in 
figure 1. The first and the last components of the system are springs that are connected to 
fixed points. 

(a) Nominal spring-mass systems 

From Newton's second law, the equations of motion of the masses are easily found to be 
mjd 2 Xj(t) 

1 d J = kjX hl (t) - (kj + k j+1 )xj(t) + k j+1 x j+1 (*) +fj(f) for ; = 1, . . . , N, (2.1) 

where kj denotes the spring constants, ray denotes the masses, and xj(t) and fj(t) are, respectively, 
the displacements of the masses and the external forces as a function of time, t. Note that 
£ 0 (f) = X]v+i(£) = 0 is defined for conciseness. The Fourier transform gives the frequency domain 
equations that may be expressed in compact matrix form as (K — co 2 M)x(co) = i(co) with diagonal 
mass matrix M, tri-diagonal stiffness matrix K, and input and output vectors i{co) and x(co), 
respectively. 

The nominal system considered here has all its masses, as well as all its spring constants, equal 
to each other, kj = k,j = 1, . . . ,N + 1 and ray = m,j = 1, . . . , N. Based on this assumption, the matrix 
K is k times the tri-diagonal matrix T, which has 2's on its diagonal and — l's on its two neighbour 
diagonals, while the matrix M is m times the identity matrix I. It is also assumed that all the 
components of the force vector have the same time dependency so that one can work with the 




constant vector, f = \f\ ... /n] T - This corresponds to an impulse force, f = f8(t), where 8(f) is the 
Dirac delta function. In matrix form, the nominal equations of motion are therefore 

(T-^(=),) W .(i)l (2.2) 

The normalized case k = m = l covers a wide range of cases. For example, the present case can 
be expressed with the normalized frequency parameter, k = a) 2 (m/k), force vector, f(k) = (l/k)f(co), 
and output vectors, x(X) = x(co). There is also no restriction to undamped systems, because the 
generally complex frequency parameter, k, permits the consideration of damping. 

(b) Disturbed spring-mass systems 

Some disturbance is allowed to the nominal spring-mass systems. Specifically, a different 
impedance is considered for the subsystem made of the two extreme (first and last) masses, as 
illustrated in figure 2. This disturbance might, for example, be an additional spring connecting 
the two masses, a small mass added on the first mass, or an additional spring connecting the 
last mass to a fixed point. Any such disturbance is a priori allowed (one does not prescribe 
the disturbance to necessarily be a combination of springs or masses or even physically possible). 
In order to conveniently study its effect on the nominal system, it is expressed as the product of 
a unit disturbance scaled by a scalar, s. In matrix form, this results in the bi-parameter system 
(T — co 2 (m/k)l — sS)x(co,s) = (l/k)f, where the three normalized matrices are 
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As in %2a, the normalized vectors are x(k, s) = x(co, s) and f = (1 /k)f for k = co 2 (m/k). For simplicity, 
S is rank-one, i.e. it may be expressed as the outer product of two vectors, S = sisj". 

(c) Normalized spring-mass systems 

It was seen in %2a,b that the equations of motion can be written in a normalized form that is 
equivalent to working with unit mass and stiffness, and with a frequency parameter, k = co 2 . The 
same normalized systems, (T — AI)x(A, 0) = f and (T — kl — sS)x(k, s) = f , are studied in the rest 
of the paper. They support a wide range of situations as illustrated in the following subsections. 
Note the mnemonic notation, T for fri-diagonal, I for identity and S for square. 




(i) First application: normalized damped system 

A first normalization that involves damping is now illustrated with the damped benchmark problem 
schematized in figure 3. The nominal equations of motion are found to be 

d 2 Xj(t) _ fdx j+1 (t) dxj{t) dx hl (t) 
= (f) ~ *W + + c {-^ 2-£- + 

d£(f) - 

-kaX,<t)-C a -± r +fj&(t), (2.4) 

where c is the additional damping coefficient and k a and c a are the real positive stiffness and 
damping coefficients of the additional local elements. This results in the nominal matrix system 



T (co 2 m-icoc a -k a ) 
(k + icoc) 



(k + icoc) 



with complex frequency argument k = (co 2 m — \coc a — k a )/(k + icoc) and normalized force 
\/{k — icoc)i. Similarly, considering the disturbance leads to a non-normalized system such that 

[(k + icoc)T -(oo 2 m- iooc a - k a )\ - s(k + icoc)S]x((o, s) = f . (2.6) 

Particular case. Equation (2.4) may appear in many contexts. It can, for example, describe the 
linearized behaviour of nonlinear systems around an equilibrium or solution point. Particular 
cases of such a situation include that of the linearized discrete form of the sine-Gordon 
equation d 2 0/d£ 2 — d 2 0/dx 2 + sin(0) = 0 that notably appears in the modelling of magnetic flux 
transmission of Josephson-junction transmission lines [26]. The poles of this linearized sine- 
Gordon equation around soliton solutions could then provide information about the spectrum 
of the phonons and the corresponding eigenmodes would be the phonons far away from the 
solitons [27,28]. Considering equation (2.4) and disturbed spring-mass systems, such analysis 
could possibly involve forces, damping and disturbance. 



(ii) Second application: normalized multi-degree-of-freedom system 

The second normalization example and benchmark application involve non-normalized springs 
that are multi-degree-of-freedom (multi-d.f.) components. The particular model considered is 
schematized in figure 4. Each consecutive pair of the N masses is separated by an equivalent 
or generalized spring that consists of two masses, m\ and mi, and four springs, k\-k\. The full 
system has therefore more than N degrees of freedom. The transformation into the normalized 
form starts with the resolution of the equations of motion inside individual spring components. 
Looking specifically at the component between the (/)th and (j + l)th masses and at mass j gives 

- kiXj_i{oo) + (ki +k 2 + h- oo 2 mi)yf{oo) - k 2 yf(oo) - k±Xj{oo) = 0, (2.7) 

- k 2 yf(co) + (k 2 +h- co 2 m 2 )yf{co) - k 3 Xj(co) = 0 (2.8) 
and (h +h + U-oj 2 m)Xj(co) =fi + hyf(oo) + k 3 yf(aj) + hy^ (co). (2.9) 



Figure 4. Illustration of a regular collinear system with multi-d.f. springs. 



Using equations (2.8) and (2.7), the internal variables, yf{co) and y^{(o), can be solved in terms of 
the neighbouring variables, Xj_\ and Xj, and substituted in equation (2.9). Basic algebra gives, for 
identical mx—m^ — m s and k\ = ■ • • = k^ = k s , and with co s 2 = k s /m s , 



(/')/ x [(2-co 2 /a) s 2 )x hl (oj)-(3-oj 2 /oj s 2 )]xj(a)) 
[1 - (2 - a>*/cD£)(3 - o) z /o)i)] 

j ,(/) ^ ^ _ _ [g-lH - (4 - ft) 2 /ft) s 2 )*;M] 

[1 _ (2 - a; 2 /^s 2 )(3 - ^ 2 )] 
.,2 2x 



yi\o)) = - J \ ' (2.11) 

2 - /,.,4V 



and f3 - co 2 ^\ x ; H = £ - - (3 ?/ ^ s2) 7/ 2xi [*/-iM + 3*/M + ^lHl- (2-12) 

V k s J 1 k s [1-(2-oj 2 /oj 2 )(3-co 2 /oj 2 )] 1 1 1 

Although it may not be obvious, this equation can easily be written up in the normalized form 
by two operations: the multiplication of both sides by u(co) = — [1 — (2 — co 2 /o) 2 ){3 — co 2 /oj s 2 )]/ 
(3 — co 2 /co s 2 ), and the addition and substraction of 5xj(co) from the bracket [xj_\ + 3xj + xy + i] = 
[xj_i — 2xj + Xj + \\ + 5xj. This results in the desired form (T — AI)x(A) = iu(X), where X = 5 — u(co) 
(3 — o) 2 m/k s ) is a rational function of a* 2 with cubic numerator and linear denominator. For each 
value of A equal to a pole or zero of the normalized system, it is thus possible to retrieve the 
corresponding values of co 2 by trivially solving a cubic polynomial. 

Some care must now be taken with the right-hand side, i.e. with the normalized vector, 
fu(X). Its scalar magnitude can indeed be zero or infinite for discrete values of co 2 that cancel 
the quadratic numerator or linear denominator of u(X). Possible simultaneous cancelling of the 
numerator and denominator must be carefully considered as well as the a priori possible cases 
where a pole or zero of the system (T — AC) is neutralized, respectively, by a zero or pole of the 
magnitude, u(X). With these considerations, the zeros and roots of the normalized system and of 
the scalar function u provide all those of the non-normalized system. Note that in the case of the 
first benchmark, i.e. in equation (2.5), the frequency-dependent scaling function l/(k + icoc) can 
be neither zero nor infinite as long as the real k and c are strictly positive. All zeros and poles of 
the non-normalized and normalized systems (for a constant force) are therefore identical. 

For the current multi-d.f. benchmark, special attention must also be paid to the additional 
internal d.f. In order to avoid missing internal modes, the poles of the internal systems — in the 
present case, the roots of the polynomial [1 — (2 — co 2 /co 2 )(3 — co 2 /co 2 )] — must be analysed. Finally, 
the study of the transfer functions and zeros of the normalized system provides direct information 
only for the main masses of the whole original system. More information is however available. 
For example, equations such as (2.10) and (2.11) provide all the system responses from those at 
the main masses. These expressions can thus be used to carry out the work of identification of 
zeros at the internal masses. By reciprocity, the unit transfer functions from an internal mass to a 
main mass can also be treated in a dual manner. 

Other cases. Under the same umbrella as that of the current benchmark, multi-d.f. systems such 
as in figure 4 allow one to deal directly with other systems such as the one illustrated in 
figure 5, where a modified stiffness component (the 1 d.f. subsystem circled by a dashed line) 
only appears every other spring. In this case, this is simply done by using k\—k and m\ — m. 
The case of diatomic or multi-atomic lattice vibrations [2] is also covered by this benchmark, 
for both alternating masses and alternating springs. Periodic composite materials such as those 



Figure 5. Illustration of a colinear spring-mass systems with alternating regular and one-degree-of-freedom springs. 



described in Silva [29] and Andrianov et al. [30] can be put in the same framework. Furthermore, 
the present approach can be directly applied to condense macro cells of a homogeneous chain, 
when considering a so-called multi-field approach of vibrations as in Vasiliev et al. [31]. This can 
be illustrated by the case of figure 4 with m\—m2 — m, k\—k2 — k^ = k and k± = 0. 

(iii) Third application: bubbles vibrating in an acoustic field 

The third application concerns the case of gas bubbles in water vibrating in an acoustic field. The 
dynamics of such bubbles is a very complex and important field of study, as its behaviour is highly 
nonlinear for large radial changes [32] so that it can lead to locally extreme conditions of velocity, 
pressure and temperature. The study of the linear vibrations of the bubble surfaces is also very 
important as their acoustic response, when submitted to an acoustic field, can be used as a tool to 
identify and count individual bubbles from a gas leakage [33]. In the presence of several bubbles, 
there is interaction and resonance because of coupling through radiated acoustic pressure from 
one bubble to the other. This effect can be modelled by this theory for the bubble vibrations in 
their linear regime. 

The corresponding normalization is demonstrated here by starting from the coupled linearized 
equations of motion of two identical bubbles [34-36] 

^(f) + + hi ^(f) + te^jrW + co 2 0 n{t) = p ac (f) (2.13) 

and 

§ W + + h2 ld {t) + 8ll ld {t) + w l r ^=V*S), (2.14) 

where Tj denotes the variation of the radius Rj of the ;th bubble compared with its equilibrium 
value, Ro, i.e. Rj(t) = Ro + Tj(t); a = Ro/d is the ratio of this equilibrium and the distance, d, 
between the two bubbles; p ac is a force acting on the system that is proportional to the incident 
acoustic field in which the bubbles are residing; and coo is the natural Minnaert's frequency of 
resonance of the bubbles. The factors 8ft are coupled damping coefficients that are here chosen by 
symmetry, such that 812 = 821 = 8 C and 822 — hi = <$a- 

The first step to normalization is the premultiplication of the system, expressed in matrix form, 
by the inverse of the matrix a = [ * ^ ] . In the frequency domain, this gives 

( K T-a> 2 n\)r=^[l 1] T , (2.15) 
1 + a 

where the frequency-dependent equivalent stiffness and mass of the bubbles are k = co^ 
(a/(l - a 2 )) - - a 2 ))(8 12 - a8 22 ) and fi = (1/co 2 )[2k - (oo 2 /(l - a 2 )) - (ia>/(l - a 2 ))(8 n - 

01821) + co 2 ]. The division of both sides of this equation by k then provides the normalized equation 
with X = co 2 /jl/k and right-hand side f = [1 l] T /? ac /[/c(l + a)]. 

All analysis in this paper can therefore be directly applied to analyse the vibration of gas 
bubbles in liquid. In particular, all the results about the transfer functions, guv, from the first to the 
last mass correspond to the transfer functions between the two bubbles in the present case, N = 2. 
Note that even though the boundary conditions of the bubbles might not formally be equivalent 
rigid attachments, this model can perfectly be used, without any approximation or error, to 
analyse the two bubbles. One advantage of the theoretical model of this paper is that otherwise 



(a) (c) (d) 




Figure 6. Schematic of some fields of application of this work: (a) chain of molecules in a single-walled nano-tube; (b) vibrating 
diatomic lattice with alternating springs or alternating masses; (cj vibration of a bladed disc and its simplified model as a circular 
vibrating system; and (d) acoustically excited gas bubble in liquid. 



difficult to measure properties of the bubble-fluid system can be evaluated: by matching the 
predicted loci of the resonances and zeros to their evaluations obtained in a well-controlled 
experimental set-up, the actual values of parameters such as damping can be appropriately tuned. 
In the context of identification of gas leaks from pipes, natural methane seeps, or from undersea 
carbon capture and storage facilities [33], one can also push the analysis further to assess how the 
cross-bubble coupling affects the acoustic scattering from an ensemble of bubbles and how the 
consideration or non-consideration of this coupling impacts the measurements of the number of 
bubbles and quantity of leaking gas. The vibrational behaviour of pairs of bubbles with regard 
to that of single bubbles is also important in shock wave lithotripsy [37] as the source image of a 
bubble close to the surface of a kidney stone is exposed to a symmetric source image. 

General applicability. As illustrated in the previous sections and in figure 6, the work presented 
here can generally be applied and extended to various fields and the d.f. are not restricted to 
displacements. They may, for example, be rotations of a beam or values of the acoustic pressure at 
different points of a waveguide. Note that figure 6c corresponds to a circular chain with damage 
whose transfer functions have been studied in [13, section 5, appendix D], starting from this 
theory. 

3. Exact expressions of the nominal and disturbed transfer functions 

Considered in this section are the responses or transfer functions between a unit force applied at 
the kth mass and the displacement of the ;th mass. Such transfer functions are denoted 

g j/k (Ks) = ej(T -XI- sSrV, (3.1) 

where ey is a unit vector with only non-zero coefficient at position j, in the general, disturbed, case 
while the nominal transfer functions in the particular case of absence of disturbance are 

gjik (X,0) = ej(T-Xl)- 1 e k . (3.2) 

The transfer functions for general input and output vectors, f and c, are available from linear 
combinations of those as c T (T — kl — sS) -1 f = J2j=i, ,N J2k=i, ,N c j8j,k(^' s )fk- Both the nominal 
and disturbed transfer functions exhibit remarkable characteristics that are presented further in 
the paper. Exact expressions of the transfer functions are first derived in §3b,c, by making use of 
their link with Chebyshev polynomials which is first explained in §3a. 

(a) Transfer functions and Chebyshev recurrence relations 

The nominal transfer functions have been studied at least since the time of Rayleigh. They have 
interesting properties by themselves because of the structure of T, which, for a fixed value of 



N and k, and even with (fixed) non-zero disturbance, has the inherent recurrence relation of 
Chebyshev polynomials, presented here in the general case, 

gj+l,k(K s) = (2 - X)g j/k (k, s) - g hU (K s) for ; # 1, k, N. (3.3) 

Such recurrence relations are standard in the generation of orthogonal polynomials (e.g. [38]). 
They are in particular the recurrence relations for the Chebyshev polynomials of second kind, 
U n (x), which are defined by 

U j+1 (x) = 2xUj(x) - U hl (x) for; = 1,2, . . . (3.4) 

together with the initial functions, Uo(x) = 1 and U\(x) = 2x. An alternative trigonometric form of 
these polynomials is [38, eqns (22.3.16)] 

^(cos(,))=^±^. (3.5) 

Note that equation (3.3) differs only from equation (3.4) by a shifting and scaling of the 
argument (x = l — X/2). This is the same kind of operation as the shift used to obtain 
the shifted Chebyshev polynomials (e.g. [39, p. 179]). The recurrence relations and the definitions 
of Chebyshev polynomials allow easy derivation of the transfer functions, as presented in the 
next sections. 

(b) Explicit expression of the nominal transfer functions 

Explicit expressions of nominal transfer functions are summarized in the following lemma. 

Lemma 3.1. Consider the tri-diagonal dynamic system 

(T-kI)x(k) = e k , (3.6) 

where T is the tri-diagonal matrix of dimension N with all twos on its diagonal and all negative ones on its 
neighbour diagonals, I is the identity matrix of dimension N and e k is the kth unit vector. 
Then, the components, gj /k (k, 0), of the response, x(k) = [g\ f k(X, 0) . . . gN,k(k, 0)] T , are 

^ (X '° )= U N (l-X,2) ' t0r ' = 1 k (3 - 7) 

and 

U N _ ; (1 - V2)l4-l(l - X/2) . 

^ {K0) = U N (l-X/2) ' f ° r ' = fc (3 ' 8) 

where the function Uj(x) is the Chebyshev polynomial of second kind of order j (Ho 00 = 

1, Ui(x) = 2x, Uj+\(x) = 2xllj(x) — Uj_i(x), for j = 1, . . .). By symmetry, gj^iK 0) = gk,j(K 0) m & 
gN+i-j,N+i-k(K 0) = gj,k(K fyfor all X, N,j, k. 

Proof For s = 0 and fixed values of N and k, one uses the notation Xj(k)=gj f j c (k / 0) for 
compactness. The recurrence relations above and below row k of equation (3.6) give 



xy+i(A.) =2(l-^\ Xj(k) - x hl (X) for ; = 2, . . . ,k - 1 



(3.9) 



and 

x hl (X) = 2(l-^j Xj(k) - x j+ i{\) for; = k + 1, . . . ,N - 1. (3.10) 

Considering first 2 < k < N — 1, the first and last rows of equation (3.6) may be written 
*2(A.) = 2(1 — X/2)x\(X) and x^-\(X) = 2(1 — X/2)x^(X). Along the recurrence relations, they give 

x j (k) = x 1 (k)U hl (l - 0 forj = l,...,k (3.11) 

and 

Xj(k) = x N (k)U N -j (^-^j tor j = k,...,N (3.12) 



Table 1. Samples of the nominal transfer functions from lemma 3.1 for the normalized system. 



Dim. 
N 


input 

k 


output 

J 


transfer function 






poles 


zeros 


1 
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1/(-A + 2) 






2 




2 


1 


1 


(-A + 2)/(A z - 


4a + 3) 




1,3 


2 




1 


2 


1/(A 2 - 4A + 3) 






13 




3 


1 


1 


(X 2 - 4A + 3)/(- 


-A. 3 + 6A 2 - 


10A. + 4) 


0.5858,2,3.4142 


13 




1 


2 


(-A+2)/(-A 3 


+ 6A. 2 - 10A 


+ 4) 


0.5858,2,3.4142 


2 




1 


3 


1/(-A 3 + 6A 2 - 


10A + 4) 




0.5858,2,3.4142 






2 


2 


(A 2 - 4A + 4)/(- 


-A. 3 + 6A. 2 - 


10A + 4) 


0.5858,2,3.4142 


2,2 



and, in particular, the expressions of x\ and xjsf in terms of x\, x\{X) = Xfc(A)[l4_i(l — A/2)] -1 
and X]v(A) = x^(A)[LJ]v_fc(l — A/2)] -1 . Substituting these in equations (3.11) and (3.12), at 
j = k—l and j = k + l, respectively, and carrying the results in the kth row of equation (3.6), gives 



x k (k) ■- 



[- 



U t _ 2 (l - A/2) 
U fc _l(l-X/2) 



+ 2 1- 



U N -t-i(l-A./2) ' 
U N _ fc (l-X/2) J 



(3.13) 



Grouping the two first terms in brackets on the right-hand side of this equation and using the 
recurrence relation of Chebyshev polynomials results in 



x k (X) ■- 



14(1 -A/2) !%_*_!(!- A/2) 



U fc _i(l - A/2) U N . k (l - A/2) 



(3.14) 



Combining the two terms in brackets on the right-hand side of this equation, considering the 
numerator, and using the following relation k — \ times (with q = k,k— 1,...,2 and x = 1 — A/2): 



U q (x)U N - q (x) - U q -i{x)U N - q -i{x) = U q -i{x)U N - q+ i{x) - U q - 2 {x)U N - q {x) f 
and then the recurrence relation of Chebyshev polynomials a last time, gives 

lZ N _fc(l-A/2)l4_!(l-A/2) 



*fc(A) = 



IM1 - V2) 



(3.15) 



(3.16) 



which can then be combined with equations (3.11) and (3.12) to give the expected results for 
k = 2, . . . ,N — 1. Note that equation (3.15) can be obtained by solving for 2x in the two occurrences 
of equation (3.4) obtained for; = N — q and q — 1, respectively 

For k = 1, the first row of equation (3.6) and equation (3.12) (which also holds for k = 1) gives 
xjv(A) = l/(U]v(l — A/2)) and the expected result. The same approach holds for k = N. The two 
symmetry properties can be readily verified from equations (3.7) and (3.8). ■ 

It will be shown in §4 not only that exact expressions of the poles and zeros of the transfer 
functions can be evaluated from the expressions of the transfer functions but also that their 
properties can advantageously be studied directly 

Explicit rational expressions of the transfer functions obtained by making use of lemma 3.1 
are provided in table 1. While poles and zeros can be evaluated from the expressions of 
the numerator and denominator polynomials, this becomes unfeasible or unpractical for large 
dimension systems. Explicit expressions of the poles and zeros will be provided in %^a,b and will 
clarify some of their properties. A particular case of the transfer functions from first to last masses 
that will be further studied is highlighted in the following corollary 



Corollary 3.2. Equations (3.7) and (3.8) collapse to gk,i(K0) = gi,fc(A. / 0) = U^_i/U^ r gk,N(Kfy = 
gN,k(K0) = U k _ 1 /U N/ gi,i(K0) = gN,N( k >Q) = u n-i/Un, and gN,i(K0)=gi,N(K0) = l/U N where, 
for brevity, Uj denotes the value Uj(l — k/2). 



(c) Explicit expression of the disturbed transfer functions 

The following theorem provides the disturbed transfer functions for any values of s, N and k. 

Theorem 3.3. Consider the dynamic system (T — XI — sS)g^(A,s) = e^, where T is the tri-diagonal 
matrix of dimension N with all twos on its diagonal and all negative ones on its neighbour diagonals, I is the 
identity matrix of dimension N,S = s\s^ is a rank one disturbance matrix, such that s r = s r iei + s^e^ 
and s\ = snei + s^eiv, and ey denotes thejth unit vector. Denote a nm = s\ m s rn for n,m = 1,N. 
The components, gjjdK s )/ of the response vector, gk(K s ) — \gi,k(K s ) • • • gN,k(^r S )] T / are 

1 

gj,k(Ks) = - {U N _ k Uj-l + s[a n (U N _kU N _j - U N _ k Uj_ 1 U N ^i) 

+ ctNN(Uj-iUk-l ~ U N _ k Uj_iU N -i)ai N (U N _jUk-i - U N _ k Uj_ 1 )]}, (3.17) 
forj = 1, . . . ,k, and 

gj,k(Ks) = * — {Ujsf-jUk-i + s[an(U N -jU N -k - U N -jUk-iUN-i) 

+ <XNN(Uj-iUk-i ~ UN-jUk-iUN-i) + <xm(Uj-iUN-k ~ UN-jUk-i)]}, (3.18) 

for j = k,...,N, and, in both cases, Q(A,,s) = — s[(«n + oinn)Un-i + ( a iN + «ni)]- 

Here, for brevity, Uj denotes the value Uj(l — k/2) where Uj(x) is the Chebyshev polynomial of second 
kind of order j (Uo(x) = 1, Ui(x) = 2x, LT ;+ i(x) = 2xllj(x) - Uj-i(x),jbrj = 1,.. .). 

Proof. The disturbed transfer function may be expressed as a rank-one update of the nominal 
transfer function (details of this may be found in [13,25]), which gives 



gj,k(Ks)=gj /k (KO) + s 



(Sngp. 0) + S]NgjM k > Q ))( s Tlgl,k(K 0) + SrNgN,k(K 0)) 

1 - s[angi / i(A.,0) + a 1N g 1/N (X,0) + a N1 g N/ i (A, 0) + a NN g N/N (k,0)] 



(3.19) 



The nominal transfer functions are available from lemma 3.1. For the case 1 <j <k <N, one 
has gj f k(K0) = Ujsf-kUj_i/U^j and the relations of corollary 3.2. Substituting these expressions in 
equation (3.19) gives equation (3.17). Equation (3.18) is obtained similarly. ■ 

Particular cases appear. If j < k, one sees from expression (3.17) that the coefficient of an 
disappears if ; = 1, while those of cx^n or a lN disappear if k = N or j = k, respectively. Similarly, 
for ]>k,a\\, aim or a iN disappear if, respectively, k = 1, ; = N or ; = k. 



4. Properties of the nominal transfer functions 

The expressions of lemma 3.1 are completely general for any regular collinear system that has 
identical point elements and connectors between two neighbour elements, independently of their 
physical nature or the real or complex value of their parameters. 

The availability of the exact expression of the responses, particularly in terms of Chebyshev 
polynomials, allows the identification of properties of the nominal system. Some of these are now 
presented and then illustrated on the damped benchmark and multi-d.f. problem introduced in 
§2c(i),(ii). A particular focus is on the position of poles and zeros. 



Ql I I I I ULJ i i I i i i i i i 

-1.5 -1.0 -0.5 0 0.5 1.0 1.5 -1.5 -1.0 -0.5 0 0.5 1.0 1.5 
real part of frequency, real (&)•) real part of frequency, real (co) 

Figure 7. Loci of the {a) five pairs of roots and [b] three pairs of zeros of the damped benchmark for N = 5, k = 0.7, m = 1.3, 
c = 0.4, k a = 0.5 and c a varying from 0 to 2.15. The circles, crosses and squares indicate the locations, respectively, for c a = 0 
and the two smallest critical values of the poles, c a = 1-7837 and 2.0980, while the stars indicate the location of the zeros at 
their first critical value, also at c a = 2.0980. 



(a) Location of nominal poles 

For the poles of the nominal system, one has the following result, which can be, and has been, 
derived by various methods. Notably, the poles are the eigenvalues of the tri-diagonal matrix T 
and they are the roots of a shifted Chebyshev polynomial [38, eqn 22.16.5]. 

Corollary 4.1. Using the same notations as in lemma 3.1, for a given value of N, the poles of 
any transfer function, g(X,0) = c T (T — X\)~ l i, share the same poles which are the roots of the shifted 
Chebyshev polynomial of second kind and order N, 11^(1 — X/2). In particular, these roots, Xj = 2 + 
2 cos (jit /(N + 1)), j = 1, . . . , N, are real and all located within the finite range 0 < X < 4. 

Proof Trivially, from lemma 3.1, since all g(X,0) are rational polynomials with denominator 
[L7]v(l — A/2)]. The roots, cos (Jtt /(N + 1)), of U]v(.) can easily be obtained from equation (3.5). ■ 

The fact that all poles are known and real can be used to derive properties of any particular 
system that may be expressed as the nominal system. This is illustrated on two applications. 

(i) Damped benchmark 

In the case of the damped benchmark problem, the Chebyshev polynomials have argument 

l _ X = (k + k a /2) + Mc + c a /2) - co 1 m/2 
2 k + icoc 



[(k + k a /2 - co 2 m/2)k + co 2 (c + c a /2)c + ico/2[c a k - k a c + co 2 mc]] 



(4.2) 



(k 2 + co 2 c 2 ) 

Inverting this relation, the characteristic equation o) 2 m — i(Xc + c a )co — (Xk + k a ) = 0 follows. There 
are therefore two possible values of co corresponding to any value of X, 



i(Xc + c a ) ± ^(Afc + k a )m - (Xc + c a ) 2 



(4.3) 



As all poles, in the X domain, are in the open real interval 0 < X < 4 and since (Xjk + k a )m> 0, their 
corresponding values (o\ t 2 (Xj) in the co domain have a strictly positive imaginary component as 
long as c or c a > 0. The critical damping, c or c a , can be identified by the value (Xc + c a ) cr it = 
2y/[(Xk + k a )m\ that separates regions where the argument of the square root is positive or 
negative. When the argument is negative, the co-eigenvalue is overdamped, i.e. purely imaginary 
so that the corresponding modal component of an impulse response is non-oscillatory in nature. 
In general, the critical damping values vary from one eigenvalue to another. This is illustrated by 
the root loci of figure 7 for N = 5 and varying values of the additional viscous parameter, c a . 



The poles for c a = 0 are co 1/2 (kj) = {±0.7261 + 0.0412i, ±0.9484 + 0.1538i, ±1.1691 + 0.3077i, 
±1.3368 + 0.4615i, ±1.4368 + 0.5742i}, and the corresponding critical values are c a = {1.7837, 
2.0980, 2.3432, 2.4770, 2.5302}. 

(ii) Multi-degree-of-freedom benchmark 

Turning to the particular multi-d.f . benchmark, there is a total of 3N + 2 masses for N 
main masses. For N = 3 and 11 masses, there are three roots of the normalized system, 
k = 2 + 2 cos (jjt/(N + 1)) for ; = 1, . . . , N, i.e. A = 0.5858, 2.0000, 3.4142. Then, solving k = 5 - 
w(&>)(3 — a) 2 m/k s ) for a; 2 provides nine root values of oj 2 . The consideration of the additional roots 
of the polynomial [1 — (2 — co 2 /co 2 )(3 — oj 2 /o) 2 )], as discussed in %3b, provides a supplement of 
two poles, co 2 = 1.3820a; 2 , 3.6180a; 2 , that correspond to k = 5. The last possibility of resonance 
cancels the denominator of the right-hand side of equation (T — kl)x = iu at co 2 = 3co s 2 . For 
particular values of m = 2, m s = 3 and k s = 0.7, this corresponds to a degenerate equation 
with k = oo. There is no additional resonance as expected, and the whole set of 11 poles is 
oj 2 = {0.0313, 0.1126, 0.2039, 0.3225, 0.7381, 0.7496, 0.7709, 0.8442, 1.2419, 1.3545, 1.4472}. 

(b) Location of nominal zeros 

Contrarily to the poles that are identical for any transfer function, the location of zeros varies 
from one transfer function to the other depending on the particular input and output vectors. The 
number of zeros itself may vary from as much as N — 1 to as little as zero. Using lemma 3.1, some 
properties of the zeros of the mass to mass transfer functions, gydK 0), are now derived. 

Corollary 4.2. Using the same notations as in lemma 3.1, for a given value of N, the zeros of 
the transfer functions, gjjdKQ) = ej(T — AI) -1 e^, are the union of the roots of two shifted Chebyshev 
polynomials of second kind, U mi (l — k/2) and U mN (l — k/2), of respective orders m\ =min(/ — 1, 
k — 1) and m^ = min(N — j,N — k). These zeros, {k^ = 2 + 2 cos (mn /(mi + 1)), m = 1, . . . , m\ } and 
{k^ N ^ = 2 + 2 cos (mn /(m^ + 1)), m = 1, . . . , m^}, are real and located within the finite range 0 < k < 4. 

Proof. Trivially from the lemma and the evaluation of the roots of the two Chebyshev 
polynomials at the numerator of the transfer functions. ■ 

The existence and location of the zeros appear to be slightly more elusive than those of the 
poles. This is particularly true when general input and output vectors, f and c, are selected. Even 
if such vectors were real, there would be no insurance that the zeros of the transfer function 
g(k,0) = c T (T - AJ) -1 f would be real. This might seem a contradiction as this transfer function 
would be a linear combination of real polynomials with real poles and zeros, but is simply 
explained by the fact that any given real polynomial of order N — 1 or less may be generated 
as a real linear combination of the numerator of the nominal transfer functions ^(A,0) = ej 
(T — AI) -1 efc, and that polynomials with real coefficients may have complex roots. 

The general case does not preclude the existence of properties in particular cases. In 
corollary 4.2, for unit input and output vectors, one has full information on the existence and 
location of zeros. This includes the following results. 

Corollary 4.3. Corollary (4.2) shows that the zeros of the nominal transfer functions g^ have double 
multiplicity when j — 1 = N — k. Higher multiplicities are impossible. 

Corollary 4.4. Corollary (4.2) shows that there is only one finite zero of the transfer function between 
the first to the penultimate masses of the system. Furthermore, this zero is simply at frequency k = 2. 

As for the case of poles, results of corollary 4.2 can be applied to any system that may be 
expressed in the form of the nominal system. This allows one to draw conclusions on the location 
of zeros for a wide range of problems. For example, in the case of the damped benchmark 
problem, one can use the inverted relation (4.3) and the discussion that follows to prove that 
all zeros of this problem are complex, with a strictly positive imaginary component as long as c 
or c a > 0. 



Additional information may be obtained from Chebyshev polynomial properties. 

(i) Damped benchmark 

In the case of the damped benchmark for N = 5 and considering the transfer function from the 
second to the third masses, #2,3/ there are three a priori zeros, which are plotted in figure 7, as 
m\ — \ and = 2. One can, however, note that each pair of zeros coincides with and therefore 
cancels a pair of poles. This corresponds, in the case of the first zero, to the value (j = 4)/(N + 1), 
which equals (j = 2)/{m^ + 1) = 2/3. The cancelling of the poles by a zero in the transfer function 
between the second and the third masses for N = 5 is a general property that is independent of 

the parameter values, k, c, k a , The critical values at which the zeros have or not a real part are 

defined by the same equation in terms of k as that for the critical values of the poles. 

While all zeros of all transfer functions from one to another mass coincide with poles in the 
case N = 5, its other transfer functions #2,2/ #3,3/ #4,4/ #2,4/ #4,2 still exhibit actual zeros because of 
their double multiplicity. The case N = 5 is somewhat particular, as, in general, not all zeros of 
unit transfer functions coincide with poles. For example, there is no such coincidence for N = 4. 

(ii) Multi-degree-of-freedom benchmark 

Considering the multi-d.f. benchmark, again for N = 3 and a total of 11 masses, all the zeros of 
the unit transfer functions between the third, sixth and ninth masses are provided by corollary 
4.2. The only possible cases where zeros of unit transfer functions between the main masses may 
exist correspond to non-zero values of m\ and m^. This can only happen when; ovk = 2, i.e. when 
the second main mass, i.e. the sixth mass, is involved, at the location k = 2 + 2 cos(7r/2) = 2. This 
zero is cancelled by a pole for g\ f 2 and #2,3/ so that the only actual zero with double multiplicity 
is only present in the transfer function g2,2 at the sixth mass. The additional possibility of zeros 
between the three main masses has to be considered in the case where the numerator of u(co) 
cancels on the right-hand side of (T — kl)x = iu. For m = 2,m s = 3 and k s = 0.7, this happens at the 
locations oj 2 = 13820oj 2 , 3.6180&>g where k = 5. These are actual zeros, as k = 5 is not a root of the 
normalized system, and the denominator of u(co) does not cancel at these zeros. 

(c) Magnitude of the nominal transfer functions 

The focus of this section is on the particular transfer function from the first to the last masses. As 
it is the inverse of a Chebyshev polynomial, as presented in corollary 3.2, its properties may be 
obtained directly from properties of these polynomials. In particular, the following well-known 
equioscillation property leads to the fact that, in the range of real parameter k e [0, 4], the transfer 
function is outside an ellipse. 

Lemma 4.5. The scaled Chebyshev polynomial of second kind, y/l— x 2 Uj(x), equioscillates 
j + 1 times between —1 and 1 in the range xe[— 1,1]: it has (j + 1) local extrema in this 
range and their values alternate between —1 and 1 for increasing values of x. The extrema 
occur at x m = cos (((/ — m + 3/2)jt)/(j + 1)), m = 1, ...,/ + 1 and the last extremum is a maximum, 

yi-^ +1 U;(%+l) = l. 

Proof In the range x e [— 1, 1], 6 = acos(x) is real and decreases from tt to 0. From equation (3.5), 
one sees that the scaled function y/l — x 2 Uj(x) is sin((/ + 1)0), whose alternating extrema —1 and 
1 are reached at (j + 1)0 = (j+ l/2)n, (j - 1/2)tt, . . . , 1/2tv. ■ 

The scaled Chebyshev polynomials of second kind are presented in figure 8 for the few first 
orders. The reader intrigued by the 'white curves' in this figure is referred to an 'asin-acos' version 
of the plots where the curves are straight lines and to Ortiz & Rivlin [40], where this phenomenon 
is explained. Such properties of the scaled Chebyshev polynomials translate into properties of 
the transfer functions. Notably, lemma 4.5 allows one to identify properties of the magnitude of 
the giflik, 0) transfer function as shown in the below corollary of the lemma and corollary 3.2. 
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Figure 8. Illustration of scaled Chebyshev polynomials Vl — x 2 Uj(x) for/ = 0, . . . , 4 in [a), and/ = 0, ... , 29 in (b), as a 
function of x = 1 — A /2. The points corresponding to the loci of the last roots are indicated by squares in {a). 
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Figure 9. Illustration of the exclusion zoneforthe values of the nominal transfer function^ for N = 7. No values (A, 2g^[k)) 
can lie within the exclusion circle and g Ni -\ (A) is tangent to this zone at eight points marked by squares. 



Corollary 4.6. Using the same notations as in lemma 3.1, for any real value A e [0,4], the absolute 
value of the transfer function, £n,i(A,0) = e^(T — from the first to the last degree of freedom, 

is larger than yfl — (A/2 — l) 2 . In other words, in the real plane, all points (A,2g]v / i(A,0)) are outside a 
circle of radius 2 centred at (2,0). The transfer function is tangent to this exclusion region at the points 
i m = 2 - 2 cos (((N - m + 3/2)n)/(N + 1)), m = 1, . . . , N + 1. Also, this transfer function has no 
finite zero. 

Although lemma 4.5 is trivial and has been previously noted, as incidentally in [41], 
property 4.6 of the transfer function appears to not have been identified before. The exclusion 
zone of the magnitude of the transfer function g^i is illustrated in figure 9 for N = 7. For any other 
system dimension, N, there will be N + 1 tangency points between the same exclusion zone and 
the transfer function. They are to be put in the light of those of the scaled Chebyshev polynomial 
to the constant lines 1 and —1 illustrated in figure 8. Their location at equiangles of the exclusion 
zone in the A-2gN,l plane is evident from their expression in theorem 4.6. 

5. Properties of the disturbed transfer functions 

The properties of the disturbed transfer functions may be derived from their exact expressions 
given in theorem 3.3. As previously, the focus is now on the location of poles and zeros. 

(a) Location of disturbed poles 

The locations of poles of the disturbed transfer functions vary depending on the particular 
disturbance vectors d/ and d r but, for any choice of these and a given magnitude, s, they are 
identical for all transfer functions. They are characterized by a necessary zero denominator 



Q(A,s), Q(X,s) = U N (1 - A/2) - s[(a n + a NN )U N ^(l - */2) + («in + am)] = 0 and are therefore 
the zeros of a polynomial of order N in A. An alternative point of view is to consider Q(A,s) to 
be a function of s (at a fixed value of A). Two cases may happen: either both U^(l — A/2) and the 
coefficient of s are zero, in which case any value of s corresponds to a pole, or there is just one 
pole which is the single root of the polynomial of order 1 in s. The first case might only happen 
at a finite number of values of A as Uiv(l — A/2) and [(an + ajsnsi)Ujsr-i(l — A/2) + (aw + «ivi)] 
are polynomials in A of order N and N — 1, respectively, which therefore have only N and N — 1 
zeros each. An example of the first case is given below. In the second case, the poles are simply 
defined by the roots, as a rational function of A. This is expressed in the following theorem, whose 
proof is trivial. 

Theorem 5.1. Using the notations and assumptions of theorem 3.3, for any vectors c, f and values of 
a nm 's, all poles of the transfer function c T (T — AI — sS) _1 f are identified by either the function 

5 W= U N (l-X/2) 

[(a n + a NN )U N _i(l - A/2) + (a 1N + a N1 )] 

for any value of X, or any value of s at the discrete values of X, if any, that cancel both the numerator and 
denominator of expression (5.1). 

The following particular case, which includes additional springs or dampers connecting the 
first and last massses, is remarkable. 

Theorem 5.2. Using the notations and assumptions of theorem 3.3, if s r \ = — s r ]v or sn = —s\n then 
the poles, X, of the transfer function g(X, s) = c T (T — AI — sS) _1 f are equal to 

i, = 2(l-cos(4)) where 4= *^*^ , * = 0 fl (^) ( 5 - 2 ) 

for any value ofs and 

~m 1 sin(((N+l)/2)acos(l-A/2)) 

S M = - 777^t — i\~r>\ 7^ where tr(S) = a n + o^n (5.3) 

tr(S) sm(((N - l)/2)acos(l - A/2)) 

for any value of X. The 'floor' value fl((N — l)/2) is (N — l)/2 if it is an integer, or N/2 — 1, otherwise. 
Note that the value of s could be expressed as a ratio of Chebyshev polynomials of second kind with half 
orders 

§ 1 H(N-1)/2(1-V2) 

Proof. In either case of equation s t \ = — s r ]v or sn = —s\n, one has (an + cinn) — — ( a iN + <*ni)- 
The denominator Q(A, s) of the transfer functions with unit input and output vectors described in 
theorem 3.3 is therefore equal to 



1-^) 
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It is a polynomial of strict order N in A and therefore has exactly N roots which are the 
poles of the transfer functions. By definition of the Chebyshev polynomials of second kind, U^ 
(1 - A/2) = sin [(N + 1)9]/ sin (0) where 6 = acos(l - A/2) so that 

{sin [(N + 1)0] - s tr(S)[sin (NO) - sin (0)]} 
Q(Ks)= — . (5.6) 



As A is a function of cos(0), the poles can be identified by the unique values of 0 in the interval 
[0,7r] that cancel the right-hand side of equation (5.6). Particular attention is paid to the values 



0 = 0 and tt as they also cancel the denominator sin(0). Looking at the numerator, the sinusoidal 
terms can be developed as [38, eqns(4.3.24/35)] 

sin [(N + 1)0] = 2 sin (^y^) cos (^T^) (5 - 7) 

and 

sin (NO) - sin (0) = 2 sin cos (~^~ 6 ')' ( 5 ' 8 ) 

which shows that 

Q(X,s) = {2 cos (^<?) sin (^*) - str(S)sin (^^) J / sin W- (5-9) 

The cosine function is zero for all values 0\, k = . . . , — 1, 0, 1, Those, however, only correspond 

to distinct real zeros, i k , of Q for the fl((N + l)/2) values, k = 0, . . . ,fl((N - l)/2). This is due 
to the cyclicity of cos(0) and because the case k = N/2 for even N cancels both the numerator 
and the denominator in equation (5.9). In the latter case, I/Hospital's rule gives 

hm^^^(d/d6'){sin(6>)} 

This shows two things, when N = 2k is even: on the one hand, ^N/2=k — 4 is not a pole for all 
values of s but, on the other hand, the pair ^/2=k = 4,s = — l/tr(S) is always a pole — this is a 
particular case of equation (5.1) for any value N. 

Now, looking at the term of equation (5.9) in brackets, in the interval 0 < 0 < tt, each of its 
roots corresponds to a pole since sin(0) > 0 in the selected interval. If both sinusoidal terms in the 
brackets were zero at the same time, one would have 0 = 27ik/(N + 1) for k = 1, 2, . . . , fl(N/2) at 
the same time as sin [(N + 1)0/2] - sin [(N - 1)0/2] = 0, i.e. simultaneously to [38, eqn(4.3.35)] 



#) sin (H ; 



:0. (5.11) 



The sine term cannot cancel as its argument is 0 < 0/2 < tt/2. The cosine term cancels if its 
argument (N/2)0 = jrkN/(N + 1) is equal to tt/2 +/7T where ; = kN/(N + 1) - 1/2 is an integer. 
This is however impossible since kN/(N + 1) < k as long as N > —1 and kN/(N + l) = k — k/(N + 
1) > k — 1/2 since k < fl(N/2). Therefore, the term in brackets in equation (5.9) can only cancel in 
the interval 0 < 0 < tt when 

(s)= sin(((N+l)/2)fl) = sin(((N + l)/2)acos(l-A/2)) 
S ^ } sin (((N - l)/2)<9) sin (((N - l)/2)acos(l - A/2)) * K ' } 

The only values of 0 that cancel the sine terms and have not yet been considered are 0 = 0 and 
also 0 = tt if N is odd. These values however also cancel the denominator in equation (5.9) and the 
actual non-zero values of Q(k, s) are found to be the following by I/Hospital's rule: 

Q(X,s) = 2p^-str(S)^*) if 0 = 0 (5.13) 

and 

Q(A.,s) = -2 ( -str(S)^— 1 ] if 0 = tt and N is odd. (5.14) 



In both cases, it is generally non-zero, except at s = [(N + 1)/(N — l)]/tr(S), which also 
corresponds to equation (5.1). All poles have thus been identified. The transfer function with 
general vectors c and f has the same poles as those with unit vectors. ■ 

Several general properties of the location of the poles can be derived. These include locations 
where a pole will always be present and the opposite of regions from which poles are excluded. 
The location of multiple poles and the real or complex character of the poles are also properties 
important to know in control. A sample of such properties are proved in §5fl(i)-(iii). While these 
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Figure 10. Location of poles in the case of corollary 5.3, for N = 2, 4 and tr(S) = 2.4. The point /? 4 is the only pole common to 
all even values of/V. 

results can be directly applied in the situation of any systems that can be normalized, they can 
also be extended by following strategies similar to those used in the proofs. 

(i) Fixed location of poles for any system dimension 

The following corollary of theorem 5.2 identifies fixed pole locations for any system dimension. 
Corollary 5.3. Using the notations and assumptions of theorem 3.3, ifs x \ = —s t n or % = — s\^ then 

(i) for any even N = 2k, k = 1, 2, . . . , the pair (^N/2=k = 4,s = — l/tr(S)) corresponds to a pole; 

(ii) for any odd N = 2k + 1, k = 0, 1, . . . , the pair (X = 2, s = 0) corresponds to a pole. 

Furthermore, in either case, these poles are the only ones that are independent of the particular, respectively, 
even or odd values ofN. 

Proof That the pair (^N/2=k = 4, s = — l/tr(S)) is a pole for any even dimension N = 2k was noted 
in the proof of theorem 5.2 and it can easily be checked from equation (5.3). 

If N = 2k + 1 is odd, then corollary 4.1 shows that kj=k+l = 2 is a nominal pole, i.e. it is a pole at 
s = 0. Alternatively, two situations can be considered in the odd N case, using theorem 5.2. First, 
if N = 4 m + 1 for m = 0, 1, 2, . . ., equation (5.2) with k — 2m shows that X^ =2 m = 2 is a pole for any 
value of s, therefore also for s = 0. Second, if N = 4 m + 3 for m = 0, 1, 2, . . ., equation (5.3) shows 
that s(k = 2) = {sin[(m + 1)tt]/ sin[(m + 1/2)tt ]}/tr(S) = 0. 

The unicity of the invariant points for varying even or odd N is proved by examining different 
values of N. First, for odd N, there is only one possible pole for N = 1, i.e. k$ = 2, for any s. This is 
a different value from those for N = 3. 

Second, for even N and referring to figure 10, the poles for N = 2 are along two 
lines in the real X,s plane: 'L\ such that Xq = 1 for any s and 'L^ such that s = 
l/tr(S)(3 — A). (Note that, in order to show the latter, one can apply 2cos 2 (y) = 1 + cos(2y) 
and sin(y + l/2y) = sin(y) cos(l/2y) + cos(y) sin(l/2y) to U_i/2(x) = sin(l/2acos(x)) and U1/2OO = 
sin(3/2acos(x)). This gives U-y 2 (x) = [1- cos 2 (l/2acos(x))] 1 / 2 = (1/V2) (1 - x) 1/2 and Uy 2 (x) = 

(1 - x 2 ) 1 / 2 (l/V2) (1 + x) 1/2 + x (l/V2) (1 - x) 1/2 = (1/V2)(1 - x) l l 2 [{\ + x) + x]. Therefore, U 1/2 
(x)/ U-y 2 (x) =l + 2x and Uy 2 (l - k/2)/U-y 2 (l - X/2) = (3 - A).) Now, considering these two 
lines, one may consider the poles for N = 4. These happen along three lines: the vertical constant 
lines, at Xq = 0.3820 and X\ = 2.6180 for any s that intersect L 2 at the two points p\ = (X = 
0.3820, s = 2.6180/tr(S)) and p 2 = (X = 2.6180, s = 0.3820/tr(S)). The line L\ intersects the curve s(X) 
of equation (5.3) for N = 4 at the point p 3 = (X = l,s = s(l) = 1/(2 tr(S))) as acos(l -X/2) = tt/3 at 
X = 1. There is a single fourth intersection between the poles at N = 2 and N = 4, that is, the fixed 



point (A = 4, s = — l/tr(S)). In order to prove this, the expressions of s(A) are matched for N = 2 and 
N = 4, which gives 

1 sin ((3/2)acos(l - A/2)) 1 sin ((5/2)acos(l - A/2)) 



tr(S) sin ((l/2)acos(l - A /2)) tr(S) sin ((3/2)acos(l - A/2)) 



(5.15) 



sin f -acosfl - ^ = sin Qacos^l - ^)) sin Qacos^l - ^Yj (5.16) 

and 0 = 1 - cos ^2acos^l - = sin ^acos^l - , (5.17) 

where the trigonometric identity sm(a) sin(b) = [cos(a — b) — cos(a + b)]/2 was used to generate 
the third and the fourth identities. Note that the passage from the first to the second equation 
through multiplying both sides of the first equation by sin(^acos(l — A/2))sin (^acos(l — A/2)) 
is only valid if this term is non-zero. The fourth equation can be expressed as 
1 — cos(acos(l — A/2)) 2 = 0, i.e. A — A 2 /4 = 0. This results in two possible common poles: the 
actual fixed point, p\ = (A = 4,s = — l/tr(S)), and the possibility A = 0 that has to be discarded as 
it cancels the denominators of equation (5.15), while I/Hospital's rule shows that its right- and 
left-hand sides differ. One can easily verify that the three points p\ to p$ are not poles in the 
case N = 6. ■ 

(ii) Location of multiple poles 

The frequencies and values of the magnitude disturbance that lead to multiple eigenvalues are 
important information. In the context of theorem 5.2, the poles can only have multiplicity two 
at most. The location of the double eigenvalues as well as some of their properties can easily be 
obtained from the theorem as detailed below. 

Corollary 5.4. Using the notations and assumptions of theorem 3.3, if s r i = — s r ]v or sn = —s\n then 

(i) the poles have at most multiplicity two; 

(ii) there are fl((N + l)/2) poles, (A = i k , s = s(k k )) for k = 0, . . . ,fl((N - l)/2), with double 
multiplicity that are described by equations (5.2) and (5.3); 

(iii) the double multiplicity poles can only be along the two branches o/l/(tr(S)(l — X/2))for A in the 
open interval 0 < A < 4. 

Proof. The two first theses are trivial corollaries of theorem 5.2. In order to prove the third thesis, 
the expressions of s(A^) are developed. This gives 

h = 2 (l - cos (g+far) ^)) = 4 sin ((f +fcr ) ^ (5.18) 

and, using the identities sin(a — b) = sin(a) cos(fr) — cos(fl) sin(fr) and cos (tt/2 + kit) = 0, 

~ r v J sin(7T/2 + for) = 1 1 

1 k) tr(S) sin ((tt/2 + for)(l - 2/(N + 1))) tr(S) cos ((tt/2 + for)2/(N + 1)) ' 1 ' j 

Therefore, s(A^) = [tr(S)(l — A^/2)] -1 . Finally, equation (5.2) shows that A^ lies in the open 
interval. ■ 



Damped benchmark. 

The identification of the location of the multiple poles is illustrated in the case of the damped 
system application of §3a. Starting from the non-normalized system, equation (2.6), corollary 5.4 
can be directly applied provided that the frequency parameter, co, and the disturbance magnitude 
of interest, s f = s(k + icoc), are properly normalized. For example, if N = 5, there are only 
fl((N + l)/2) = 3 poles with double multiplicity at locations (A = Ay, s = s(A ; )) for ; = 0, . . . , 2 with 
Ay = 2(1 - cos (0j)), s(Ay) = (l/tr(S))sin (3<9y)/sin (2(9y), and (9y = tt(1 + 2j)/6. Explicitly, the multiple 



S' = 0.306649 + 0.136676i 




(b 



S' = 0.406649 + 0.236676i 



-1.5-1.0-0.5 0 0.5 1.0 1.5 
real part of frequency, real (co) 



-1.5-1.0-0.5 0 0.5 1.0 1.5 H 0 ' 2 

real part of frequency, real (co) 



Figure 11. Illustration of multiple pole location. Log 10 of the absolute value of the transfer function c T x(&>,5) of equation (2.6) 
for N = 5. All 10 poles are single in general cases such as for s f = 0.4066 + 0.2367i in (a) while there is a double pole, which is 
indicated by the added white box, at co = 0.7102 + 0.1566i for st = 0.3066 + 0.1367i in (b). The parameters are k = 0.7, 
m = 1.3, c = 0.4, k a = 0.5, c a = 0.3, arbitrary random input and output vectors, f and c, s r (1) = 1.5, s r (/V) = -1.5, 
S|(1) = 0.8,s,(/V) = -0.8, so that tr(S) = 2.4. 



poles are at the (X, s) locations (0.2679, 1.1547/tr(S)), (2, oo) and (3.7321, -1.1547/tr(S)). These are 
consistent with the identified property that s(A^) = l/[tr(S)(l — A^/2)] and it is clear that the pole 
(2, oo) is the only possible multiple pole with infinite value of s for any dimension N. 

Two simple operations provide the location of poles in the (co, s r ) plane: the solution of the 
following quadratic equation in co, Xj(k + icoc) = (co 2 m — icoc a — k a ) and the further scaling and 

shifting, s f = s(k + icoc), with the particular solutions co = cJf 2 : 



i(c a + cXj) ± J-{c a + cXj) 2 + 4 m(k a + kkj) 



(5.20) 



For the particular values k = 0.7, m = 1.3, c — 0.4, k a = 0.5 and c a = 0.3, this results in the four 
finite poles, (co, s) = (±0.7102 + 0.1566i, (0.7360 ± 0.3280i)/tr(S)), (±1.3852 + 0.6895i, -(0.4898 ± 
0.6398i)/tr(S)) and in the two poles at co = ±1.1325 + 0.4231i and infinite s. The existence of these 
double poles is illustrated in the maps of figure 11. 

(iii) Real and complex character of the poles 

In the context of theorem 5.2, if the argument X is real, the value of s(A)tr(S) is also real. This might 
not appear obvious because the arccosine of a real number may be complex. A formal proof of this 
statement is therefore given and it is incidentally shown that the poles are antisymmetric about 
k = 2 when N is odd. It might be worth noting that there is no (anti-) symmetry of the location of 
the poles when N is even. 

Lemma 5.5. Using the same notations and assumptions as in theorem 5.2, ifk is real, i.e. /fimag(A) = 0, 
the junction s(A)tr(S) is real, and 

v(N+l)/2 / r - - N (N+l)/2 



(V, _ y^2Ti) (N " 1)/2 _ ( w + y^23i) ( 

(|x|-y^T) (N+1)/2 + (-^(|,| + y^T) (N+1)/2 

'(ixi-y^T) (N - 1)/2 + ( -i)N( w+v ^T) (N - W2 



sM*(S) = 1 ( n-d/2 ; (n-d/2 in < 2 ( 5 - 21 ) 



for X > 2, (5.22) 



where x—l — X/2. Incidentally, if X is real, s(A)tr(S) is antisymmetric about X = 2 if N is odd, i.e. 
s(2 + (X - 2)) = -s(2 - (X - 2)) ifN is odd. 



Proof. An explicit expression of the complex value acos(x) in terms of the logarithm can 
be derived as follows. One denotes 6 = acos(x) = 0 r + \B\ where 0 r and Oj are real. As cos(0) = 
l/2(e l0 + e~ l0 ) is real its imaginary component sm(O r )(e~ 0i — e°<)/2 may only be zero if either 
0 r = kn for an integer k or 9j = 0. The latter corresponds to acos(x) real and values of k that 
are in the interval 0 < k < 4. It will be treated later. The former case corresponds to 6 complex, 
real regions |x| > 1 and intervals k < 0 and k > 4. As the cosine is a cyclic function, there 
are only two distinct cases, 0 r = 0 or 0 r = it, to consider, which, respectively, give x = cos(0) = 
(e - ^ + e 6i )/2 > 1 or x = cos(0) = —(e~°< + e 0i )/2 < — 1. These are quadratic equations in e 0i : one has 
either e 20i — 2xe 0i + 1 = 0 or e 20i + 2xe 0i + 1 = 0 whose roots are, respectively, e 0i = x± y/x 2 — 1 
and e 0i = — x ± y/x 2 — 1. One may only consider one of the two roots in each case as either gives 
the same ratio of sines in equation (5.12). One chooses the plus sign if x > 1, k < 0, 0 r = 0 and 
the minus sign if x < — 1, k > 4, 9 r = it. The imaginary part 9j of acos(x) is obtained by taking the 
logarithm of the expressions of e 0i . Thus, 

6 = acos (x = 1 - ^ J = iln (x + A 2 - l) if A < 0 (5.23) 

= 7T + i In (-x - y/x 2 - l) if A > 4. (5.24) 

Evaluating the sines of equation (5.12), sin[(N ± 1)0/2] = ( e W±l)e/2 _ e - i ( ]V± W 2 )/(2i) gives, 
for A<0, sin(((N±l)/2)^) = (l/2i)(e-^ ]V±1 )/ 2 ) ln ( x+ ^ ZI ) _ e ((N±l)/2)ln(x+V^T)), Similarly by 
splitting the terms of the exponents, one obtains for k > 4, 



[(N±1)0 
sm 



,i((N±l)/2)^ -((^±D/ 2 ) ln(-x- V^T) _ i( (N±l)/2)^ e ((^±D/ 2 ) ln(-x- V^l) 



2i 

(5.25) 



Equation (5.21) is obtained straightforwardly for X < 0 from theorem 5.2 and the fact that 

(which can be easily verified). Equation (5.22) is similarly 
obtained for k > 4 by additionally enumerating the four possible values, i, —1, — i, 1, of e 1 (N+ 1 > r / 2 
depending on the modulo m = N mod 4. The expressions are real since |x| > 1. 

The proof for the interior interval — 1 < x < 1, i.e. for the case 0j = 0, follows a similar path. 
One has x = cos(<9) = (e^ + e~ i6r )/2 and sin[(N ± 1)0/2] = (^±1)^/2 _ e - i ( N±1 ^/ 2 )/(2i). Solving 
the quadratic equation in e l0r gives 0 r = — i In (x ± Vx 2 - l) , of which a single sign, say negative, 
needs to be considered and, thus, 

/ r- n -\(N+l)/2 , j- -x-(N+l)/2 

_ sin(((N + l)/2)0) _ (^-^ 2 - 1 ) -(x-^x 2 -!) 
( )_si ^« N - 1 )/ 2 ^" (x-v^)^^ ( ■ } 

One arrives at equation (5.21) again by using (x + yj x 2 — 1) _1 = x — y 7 * 2 — 1 and at equation (5.22) 
by enumerating the values of (— 1)( N+1 )/ 2 for all values of the modulo N mod 4. The expressions 
are real since sin(.) is real for real argument. The antisymmetry is directly verifiable by equation 
(5.22) and considering that x = ±(k - 2) at k = 2 - (±l)(k - 2). ■ 

The fact that a function f(k) is real when k is does not automatically imply that the reciprocal 
is true. It is, however, the case for s(A)tr(S) as proved below. 

Corollary 5.6. Using the same notations and assumptions as in theorem 5.2, if str(S) is real, then all 
the poles of the transfer function g(k, s) occur at real values ofk. 



Proof. The proof proceeds by identifying N real poles for any real value of str(S). As g(k,s) 
is a rational function of k with denominator of order N for any fixed s, complex poles are then 
impossible. 



The first fl((N + l)/2) real poles are independent of s and defined by equation (5.2). The 
remaining fl(N/2) real poles for real s tr(S) are counted by studying the expression of s(X) given in 
equation (5.3). This is done by showing that s(A)tr(S) is made of fl(N/2) branches that are defined 
on separate real regions of the frequency parameter X and such that each has values that vary 
continuously from oo to — oo. 

Looking first at the function of 0, 3(2 - 2cos((9))tr(S) = sin(0(N + l)/2)/sin(0(N - l)/2), in 
the interval 0 < 0 < tt over which X = 2(1 — cos(0)) is real, one now shows that its fl(N/2 — 1) 
poles 6/H=2/7r/(N-l) / ; = l / ... / fl(N/2-l) alternate with its zeros oj- +) = 2jn/(N + 1), 

j = 1, . . . ,fl(N/2): consider a particular pole, oj~\ for some integer j e {1, . . . ,fl(N/2 — 1)}. One 

may express it as oj~^ = Ijn /{N — 1) = lf~^7t/(N + 1) so that the non-integer /~) equals f~) = 

j(N + 1)/(N - 1). Now, /H >j because (N - 1)/(N + 1) < 1 as long as N > -1. One also has that 
<j + 1 because, if it were not true, one would have j(N + 1)/(N — 1) >j + 1 and therefore 
2/ > N — 1 (note that there are no poles of s(.) if N = 1), which is impossible because; < fl(N/2 — 1). 
As ; </-) <; + 1, one has 0 ; ^ +) < 0 ; H < oj+l for all ; e {1,2, . . . ,fl(N/2 - 1)}, i.e. this proves that 
these single poles and zeros strictly alternate in the interval 0 < 0 <tt. 

Consequently, the value of the (A-)function s(A)tr(S) varies continuously from oo to — oo when it 
goes from any pole to the next in the interval 0 < X < 4 and there is a sign change at each pole. The 
sign change of the value of s(A)tr(S) at the poles arises because the poles are single. Similarly, the 
fact that the zeros are single implies that each branch has to cut the s = 0 axis only once. The signs 
at the left and right of each pole may be determined by checking that the value of the rational 
function at X = 0 is positive. As the function reaches first a zero when X increases, the function 
is negative at values of X slightly smaller than a pole. There are fl(N/2 — 1) — 1 branches of the 
function s(A)tr(S) that each vary continuously from oo to — oo between each of its poles. 

One can look again at the poles of the transfer function g(X, s) and summarize those that have 
already been identified in the range of real values 0 < X < 4. From theorem 5.2, one knows that 
all poles are on the branches defined by real s tr(S) = s(A)tr(S) and at X^ = 2(1 — cos (Ok)) where 
§ k = (tt(1 + 2k)) /(N + 1), k = 0, . . . , fl((N - l)/2) for any value of s. One already knows that any 
chosen real value str(S) corresponds to a pole if X is equal either to one of the fl((N — l)/2) + 1 
real values X k or to one of the fl(N/2 — 1) — 1 continuous branches going from one to another pole 
of s(A)tr(S). Considering the multiplicity of multiple poles, this corresponds to N — 2 poles. 

As the denominator Q(X,s) is a polynomial of order N in X for any value of s, there are 
only N poles of the transfer function, g(X,s). Its two last poles are below the smallest and 
above the largest poles of s(A)tr(S): the value of this real continuous function is — oo when it 
reaches its first pole from below. There is a first branch that spans at least all negative values 
of str(S) and all positive values smaller than or equal to its value (N + 1)/(N — 1) at 0 = 0 
(otherwise, there would be an excess of poles of g(X,s) for s tr(S) < (N + 1)/(N — 1)). Similarly, 
the last branch spans at least all positive values of str(S). The value of the first branch must 
therefore be larger than (N+ 1)/(N — 1) for all X < 0. From equation (5.1) of theorem 5.1, one 
knows that, as a ratio of polynomials with order of the numerator larger than the order of 
the denominator, the absolute value of s(X) tr(S) is infinite when X is infinite. The asymptote 
of the first (real and continuous) branch is therefore infinity, lim^-oo s(A)tr(S) = oo. Similarly, 
lim^oo s(A)tr(S) = — oo. Two additional continuous branches going from oo to — oo have been 
identified. For any value of s tr(S), N corresponding distinct (counting multiplicities) real values 
of X have therefore also been identified such that the pairs correspond to a pole and there are no 
other poles. ■ 

(b) Location of disturbed zeros 

The zeros of the disturbed transfer functions vary both for different values of the disturbance 
parameters (S and s) and for different input and output vectors. In the case of unit input and 
output vectors, the exact position of zeros may be found as the zeros of the numerator in the 
transfer function expressions of theorem 3.3. 




(a) 



to 

O 1) 

B ^ 



2.0 
1.5 
1.0 
0.5 
0 

-0.5 
-1.0 
-1.5 
-2.0 



(b) 





IT .' I 1 




exclusion region 

I 11 .\\ .w 


6 

,1 


1 - 

CD 
N 



-1 




0 12 3 4 
frequency parameter, A 



5 -1 



0 12 3 4 
frequency parameter, A 



Figure 12. Illustration of the location of zeros ofg Ni ^{X,s) and their exclusion zone. Here, A/ = 11, S|(1) = 0.3, S\(N) = 1.7, 
s r (l) = —0.5, s r (/V) = 0.2. The colour map in (b) represents the log 10 magnitude of the transfer function for various real values 
of the disturbance magnitudes. The location of the (real) zeros matches their prediction from corollary 5.8 as presented in (a) 
where one can note that, as expected, no zeros are present in the exclusion zone. Two particular transfer functions are plotted 
in (c) and (cf), respectively, for s = 1 and 1.6. The predicted zeros and poles at real values of A can be counted in these transfer 
functions: there are five zeros and 11 poles in the first case, and nine zeros and five poles in the second case. 



Theorem 5.7. Using the notations and assumptions of theorem 3.3, the zeros of the transfer function 
gj /k (X,s) cancel the polynomial 

{U N _ k Uj_i + s[an(U N _ k U N _j - LIjv-fcUy-iLTiv-i) 

+ am(Uj-iUk-i - LIiv-fc^-i^N-iViivC^iv-y^fc-i - U^-fcUy-i)]} = 0 (5.27) 

if] < K and the following one, if ] > k, 

{U]V-yl4-i + s[au(U N _jU N _ k - U N _jU k _iU N _i) 

+ a N ]si{Uj_iU k _i - U]v_yJifc_iJi]v-i) + a m {Uj_iU N _ k - U N _jU k _i)}} = 0. (5.28) 

The location of zeros depends in a more or less intricate way on the particular value of S. One 
may identify interesting properties in particular cases. 

(i) Exclusion region for zeros 

The following transfer function from the first to last masses is remarkable: 

,S U N (1 - A/2) - s[(a n + a N N)U N -i(l - A./2) + (a m + a m )] ' 

Its zeros are indeed excluded from a region centred at X = 2 for any real argument X. This 
property is illustrated in the colourmap of figure 12. The colours in each horizontal line of the 
colourmap represent the magnitude of the disturbed transfer function for a particular value of 
the disturbance magnitude, s. As illustrated in the figure, the theorems of location of zeros and 
poles allow one to predict the exact number of real poles and zeros for various values of s. 

Corollary 5.8. Using the notations and assumptions of theorem 3.3, there is no zero of the 
transfer function gN,i(^' s ) within a circle (1 /2k — l) 2 + (sam) 2 = 1 of the real (X,sa^\) plane. The 



curves of zeros are tangent to the circle at pairs of values, X m = 2 — 2cos(((N — 1/2 — m)/(N — l))7i), 
sa N1 = (-l) N ~ m sin(((N - 1/2 - m)/(N - l))n)for m = 1, . . . ,N - 1. 

Proof. As the zeros of the transfer function cancel the numerator 1 + sa]vil/iv-2(l — V2), they 
happen at values of X and sajvi such that sam =f(X), where f(X) is defined as f(X) = — [L7]v-2(1 — 
A./2)] -1 . If A is real, so is the value of the Chebyshev polynomial li]v-2(l — V2). One knows 
from lemma 4.5 that y 7 ! — x 2 Uj(x) equioscillates between —1 and 1 in the range x e [— 1, 1]. One 
therefore has that f(X)\ = |Jiiv-l(l - A./2)| — 1 > y 7 ! - (1 - A/2) 2 in the range 0 < A < 4. As /(A) is 
the value of sa^i corresponding to a pole for a value A, the first part of the proof is complete. The 
location of the tangency points is found from the location of extrema given in the equioscillation 
lemma 4.5. ■ 



(ii) Real and complex character of the zeros 

From theorem 5.7 and corollary 5.8, it is clear that not all zeros of the transfer function between 
the first and the last masses occur at real X, even if the disturbance sa^i is real. One can, however, 
derive properties of the imaginary part of X at the zeros from the exact expression (5.29). The next 
result gives exact expression of this imaginary component when the real part of X is equal to two 
and N is even. It also shows that there is no zero of this transfer function with real sam when 
real(A) = 2 and N is odd. In this case, all zeros correspond to purely imaginary values of sajvi- 
Incidentally, this shows the corollary that, if sa^i is a complex number with non-zero real and 
imaginary parts, there are no zeros of the transfer function with real(A) = 2. 

Theorem 5.9. Using the notations and assumptions of theorem 3.3, if g^i (X, s) = 0 and 
real(A) = 2, then 



(-l)^cosh^) = 2(-l) N / 2 ^fTT 

cosh((N- 1)0,-) / /TTT\ (N_1) , / J 

i( _l)(N-i)/2 msm) 2i(-l)( N +D/ 2 y^77 



sawi = 1 — — — -A- = ... r . — ... tt , for N even (5.30) 

cosh N- 1$ / ^ /TTiY _l ( ^ /~2T7\ 

for N odd (5.31) 



sinh((N - 1)0,0 / ri — 7\ (N-1) i ri — 7\ (N_1) ' 



where both x\ = (1 — X/2)/i = — l/2imag(A) and 0j = In ^— x z - + ^x 2 + 1^ are real The hyperbolic 
sinusoidal functions have their usual definition sinh(y) = — isin(iy) and cosh(y) = cos(iy). 

Proof. All values of X may be written X = 2 — 2ixj for some real Xj. Theorem 3.3 and consequent 
equation (5.29) give that the zeros occur when saj^n = — [I7iv-2(1 — V2)] -1 = — [lliv-2(i*z)] _1 - 
Sinusoidal expression (3.5) of the Chebyshev polynomials, liiv-2(cos(#)) = sin((N — l)#)/sin#, 
necessitates 6 = acos(ix z ) = 0 r + \B{, where 0 r and 0 Z - denote real numbers. As 

iXi = cos(0) = = cos(0 r ) + i sin(0 r ) (5.32) 

is imaginary, cos(# r ) has to cancel, which means that 0 r = tt/2 if 0 < 6 < tt is chosen without loss 
of generality. Substituting this value gives the quadratic equation in e 0i , e 20i + 2x z (e (9f ) — 1 = 0, 

whose only possible solution is the positive e 0i = —X[ + ^jxj + 1 so that 9j = In ^— x\ + ^Jxj + 1^. 

Substituting 0 = (tt/2) + iln (-x z - + ^xj + l) in equation (3.5) and using sin(y) = (e^ - e _i 20/(2i) 
and e 17r / 2 = i completes the proof with little algebra. ■ 

The results of theorem 5.9 are illustrated in figure 13 for a dimension varying from N = 2 to 12. 
The fact that a zero occurs at s = —l/a^i f° r N = 2 and any value of X can easily be verified in the 
expression of the theorem for m — 2. Also notable is that zeros occur only at real values of s for 
even N and at purely imaginary s for odd N. It is relatively easy to derive additional information 
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Figure 13. Illustration of the location of zeros oig N ^(X,s) when real(A) = 2, as predicted in theorem 5.9. Depending on the 
dimension /V, specifically on the modulo m = N mod 4, the zeros occur only at purely real or purely imaginary values of s. 
Here, a m = —0.85. The plots are sorted on the modulo, m = (a) 1, (b) 2, (c) 3 and (d) 0. There is no zero for N = 1 and 
s = -1/a/vi =1.1765 for /V = 2. 



about the locations where zeros are absent. For example, the following corollary, whose proof is 
trivial considering theorem 5.9, is another exclusion result. 

Corollary 5.10. Using the notations and assumptions of theorem 3.3, if the pair {X, s) corresponds to 
a zero of the disturbed transfer function ^n,i(^/ s )/ i- e - tf gN,i{^, s ) = 0 an d ifsotm is a complex number 
with non-zero real and imaginary parts, then the real part of the X parameter cannot be equal to 2. In other 
words, ifgN,i(Ks) = 0, the concomitance of the three following conditions is impossible: real(sa]vi) ^ 0, 
imag(saivi) ^ 0, and real(A) = 2. 

(c) Applications of the properties of the disturbed systems 

Besides the applications already discussed, two fields in which the current results and approaches 
about disturbed systems could particularly be used in future work are those of vibration 
localization [42,43] and mode veering [44-46]. 

At least two directions may be followed, regarding the analysis of vibration localization, 
namely the consideration either of the eigenmodes of the disturbed system or of its forced 
responses. Based on the assessment of localization criteria in these two sets of vectors, such 
localization could be matched to properties of the systems and disturbance. 

The question of mode veering can be addressed from an analytical point of view in the context 
of particular disturbances considered in this paper. More specifically, as all possible multiple poles 
have been identified in the case considered in corollary 5.4, the conditions for mode veering or 
mode crossing are readily available. Mode crossing happens when the mode loci for a particular 
system and disturbance type go through one of these multiple poles, while mode veering happens 
when these loci pass closely but not through these particular poles. 

6. Conclusion 

The results presented in this paper about the location of zeros and poles of nominal and disturbed 
chains of particles are useful at several levels. 

First, they can be applied directly to actual physical and engineering systems with repeated 
particles, bubbles and substructures, possibly after adequate normalization similar to that 
demonstrated here for damped and multi-d.f. systems and bubble dynamics. The exclusion 



results ensure that no eigenvalues, zeros or values of the transfer functions can be in given regions 
of the frequency and disturbance parameters. Such insurance may thus allow one to design 
structures that are safer and better tolerant to uncertainties and variance in their properties. The 
identification of the actual locations of zeros, eigenvalues, multiple eigenvalues and pole-zero 
cancellations of the systems has potential utility in several fields, particularly in, for example, 
those of control and propagation of uncertainty. 

Second, the various techniques used in the paper can be applied and extended to other 
situations in order to generate further theoretical results. The added or disturbed components 
can, for example, be considered at other locations than at the two last particles. Three kinds 
of disturbance may then also be distinguished: the low-rank one as discussed here, the multi- 
rank case where the disturbance is a sum of several outer products and the more general case 
where several disturbance magnitudes, s\,S2, . . ., and their associated matrices, Si,S2, are 
considered. The location of zeros of transfer functions and their properties can also be studied 
with different, possibly arbitrary, input and output vectors. 

Third, the exact expressions and properties described here can be used as a basis for 
fundamental characterization and treatment of general dynamic systems. For example, it was 
explained in the paper that systems with different alternating mass values can be treated as a 
normalized system. A question is therefore whether there exists a fundamental set of spring-mass 
systems to which general systems can be reduced. Such a task of categorization also necessitates 
the consideration of various boundary conditions for the extremities of the spring-mass systems. 

Finally, the exact results presented here provide problem cases that can be used as reference 
benchmarks for various algorithmic techniques to evaluate zeros and roots. For example, 
approaches to evaluate multiple eigenvalues can be validated by verifying that they can locate 
the double eigenvalues of the damped benchmark discussed in the paper. 
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